Quantum resonances and rectification of driven cold atoms in optical lattices 
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Classical Hamiltonian ratchets have been recently successfully realized using cold atoms in driven 
optical lattices. Here we study the current rectification of the motion of a quantum particle in 
a periodic potential exposed to an external ac field. The dc current appears due to the desym- 
metrization of Floquet eigenstates, which become transporting. Quantum dynamics enhances the 
dependence of the current on the initial phase of the driving field. By changing the laser field pa- 
rameters which control the degree of space-time asymmetry, Floquet eigenstates are tuned through 
avoided crossings. These quantum resonances induce resonant changes of the resulting current. The 
width, strength and position of these quantum resonances are tunable using control parameters of 
the experimental realization with cold atoms. 
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The ratchet effect, i.e. the possibility to rectify trans- 
port with the help of zero-mean perturbations, has been 
discussed in order to explain mechanisms of microbio- 
logical motility 0, and was applied to other situations 

as well t 2j, including quantum systems 3|. The_dissi- 

pationless limit of Hamiltonian ratchets 
the corresponding symmetry predictions 
cently successfully studied with cold Rubidium and Ce- 
sium atoms in optical lattices with a two-harmonics driv- 
ing and a tunable weak dissipation In these exper- 
iments, the mechanism of the Sisyphus cooling has 
been used in order to furnish initial conditions in form of 
an optical lattice: an ensemble of atoms localized in the 
wells of a periodic potential. In the momentum space 
this corresponds to a narrow distribution near the mo- 
mentum p = 0. This is essential for the observation of 
the rectification effect, since for broad initial distribu- 
tions the asymptotic current tends to zero. 

The dissipationless case can be readily implemented in 
atom optics by using laser beams which generate far de- 
tuned standing waves 0|- Appropriate time-dependent 
forces can be applied to the atoms by phase modulating 
the lattice beams 0. The reachable quantum regime of 
cold atoms becoming coherent matter waves calls for a 
study of a quantum Hamiltonian ratchet. Previous stud- 
ies of quantum ratchets were based on the kicked rotor 
model |ij til LUj , which is easily treated numerically, but 
posesses a broad frequency spectrum of the kick drive. 
The above mentioned experimental realization of a two- 
frequency driven classical Hamiltonian ratchet instead 
suggests to consider the corresponding quantum prob- 
lem of a particle moving in a spatially periodic potential 
under the influence of a two frequency ac force of zero 
mean and to search for quantum peculiarities of the cur- 
rent rectification. 

The Hamiltonian for a particle with position x and 
momentum p is 8] 



H(x,p,t) = ^- + (l + cos(x))-xE(t), (1) 



where E(t) is an external periodic field of zero mean, 
E(t + T) = E(t), J T E(t)dt = 0. For the classical case 
there are two symmetries which need to be broken to 
fulfill the necessary conditions for a nonzero dc current 
0. If E{t) = -E(t + T/2) is shift symmetric, then Q 
is invariant under symmetry 

S a :(x,p,t)^(-x,-p,t + T/2). (2) 

If E{t) = E{—t) is symmetric, then Q is invariant under 



Sb ■ (x,p,t) -> (x,-p,-t). 



(3) 



The phase space of system is mixed, containing both 
regular regions and a chaotic layer around p = 0. When- 
ever the system Q possesses any of the two symmetries 
S a and/or Sb .directed transport is forbidden inside the 
chaotic layer [4|]. 

The two frequency driving 

E(t, t ) = E x cos(w(i - t )) + E 2 cos(2w(i - 1 ) + 0) (4) 

ensures that for Ei,E 2 ^ S a is always violated. In 
addition Sb is violated for 9^0, ±7r. The appearance of 
a nonzero dc current J c h = lim^oo l/t J to p(t')dt' in this 
case is due to a desymmetrization of the chaotic layer 
structure (FiglQi) . It induces a desymmetrization of the 
events of directed motion to the right and left 12]. Due 
to ergodicity inside the layer, the asymptotic current 
is independent of the initial phase to, for initial condi- 
tions located inside the chaotic layer. With the specific 
choice of the driving Q it follows J c h{Q) — —J c h(—9) 
and J c h{8) = —Jch(0 + 71") 0- From perturbation the- 
ory it follows J c h ~ EfE 2 sin 9 0,0- An efficient sum 
rule allows to compute the average current J c h by proper 
integration over the chaotic layer 0|. 



The Hamiltonian Q is periodic in time with pe- 
riod T. The solutions \ip a (t)) = U(t,t )\ip a (0)) of the 
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FIG. 1: (a) Poincare section for the classical limit; (b-f) 
Husimi representations for different Floquet eigenstates for 
the Hamiltonian Q with h = 0.2 (momentum is in units of 
the recoil momentum, p r = hkn, with fcx = 1). The parame- 
ters are E\ = E 2 — 2, lo = 2, 9 = — 7r/2 and t — for (b-e), 
and Ei = 3.26, E 2 = 1, w = 3, 9 = —k/2 and to = for (f). 
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FIG. 2: (a) A part of the quasienergy spectrum as a function 
of the parameter 9. The symbols indicate the corresponding 
Floquet states shown in Fig^b-e) . (b) The mean momentum 
(p) a vs mean kinetic energy (p 2 ) a for the Floquet states (filled 
circles). The line corresponds to the cumulative momentum, 
P a (see text). The parameters are the same as in Fig^b-e). 
All the momenta are scaled in units of the recoil momentum 
p r . 



Schrodinger equation 

ih^\m) = H{t,ta)\1>(t)), (5) 

can be characterized by the eigenfunctions of {/(to + 
T, to) which satisfy the Floquet theorem: \tp a (t)) — 
e^ l ^r Lt \4> a (t)), \<j) a (t + T)) — \<j) a (t)). The quasienergies 
E a (— 7r < E a < it) and the Floquet eigenstates can be 
obtained as solutions of the eigenvalue problem of the 
Floquet operator 

U(T,t Q )\4> a (t )) =e- iE -\<t>M). (6) 

The Floquet eigenstates provide a complete orthonor- 
mal basis and the stroboscopic quantum state can be 
expressed as 

mmT,to)) ^Y, C ^) e ~ mEa \^(to)), (7) 

a. 

where the coefficients {C a } depend on to. For later 
convenience the integer a, which sorts the states 
\4> a ) such that the mean kinetic energy (p 2 ) a = 
1/T J Q ((/) a \p 2 \4 > a)dto monotonically increases. 

With the help of a gauge transformation, —> 
exp(-ia; f* E(f)dtf)\ip) 13], the solut ion of the time- 
dependent Schrodinger equation for the Hamiltonian l(T[l. 
may be written as 

|^(t)) = e~Ti Soih\p-A{t' ,t„)] 2 +(i+cosx)}dt' (g) 



with the vector potential A(t, to) = — ^j- sin(w(t — to)) — 
j 2 - sin(2w(i — to) + 0). Due to discrete translational in- 
variance and Bloch's theorem all Floquet states are char- 
acterized by a quasimomentum k with \<p a (x + 2tt)) — 
e lhK \4> a {x)) . Here we choose k — which corresponds to 
initial states where atoms equally populate all (or many) 
wells of the spatial potential. The wave function is ex- 
panded in the plane wave eigenbasis of the momentum 
operator p, \n) = -j=e mx . The Floquet propagator 
U(T, to) is obtained by solving the Schrodinger equation 
over a single period T for a sufficiently large set of plane 
waves \n) with n — 0, ±1, ±2, ±iV, where 2N + 1 is 
the total number of basis states taken into account. The 
numerical calculations follow the integration method de- 
scribed in Ref.0|j results are shown for N = 60, which 
do not depend upon further increase of N. 

If the Hamiltonian is invariant under 
the shift symmetry S a J2J, then the Flo- 
quet matrix has the property U(T,to) = 
[[/* (T/4, to) U (T/4, t )] T U* (T/4, t ) U (T/4, t ) 

. Here performs a transposition along the 

codiagonal of U. In our case S a is always violated. 

If the Hamiltonian is invariant under the time rever- 
sal symmetry Sb ©j then the Floquet matrix has the 
property U(T, t ) = U(T, i )* 0. That symmetry will 
be recovered for 8 — 0, ir. Then the Floquet matrix has 
an irreducible representation using even and odd basis 
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states \n) s>a = (|n) ± | - n))/ y/2 - S nfi . 

The mean momentum expectation value J(to) = 
lim^oo 1/t J to (i(j(t,to)\p\ip(t,to)) measures the asymp- 
totic current. Expanding the wave function over the Flo- 
quet states the current becomes 



J(to) =5>) Q |C^ )| 2 



(9) 



where (p) a is the mean momentum of the Floquet state 

\<Pa}- 

For the symmetric case 9 — 0, ±7r it follows that 
(p)a — for all a. Consequently J = in this case. 
We especially note that this is true for states with arbi- 
trarily large kinetic energy, for which the corresponding 
quasienergies become almost pairwise degenerated. For 
6^0, ±7r the Floquet states become asymmetric. The 
quasidegeneracies are removed (FigEk)- Especially Flo- 
quet states with large kinetic energies acquire large mean 
momenta (FiglJt) , thus becoming transporting. This is a 
consequence of the fact that the true perturbation param- 
eter regulating the desymmetrization around the sym- 
metric quasidegeneracy points for small 9 is proportional 
to nE\E 2 9. FigEt also shows the cumulative average 
momentum, P a+2 = P a + (p) a +i + {p) a +2, Pq = (p)o- The 
asymmetry stems mainly from Floquet states located in 
the chaotic layer region of the classical phase space (see 
results below). With increasing (p 2 ) a , P a goes to zero in 
full accordance with the fact that total current over the 
whole momentum space should be zero Q. 

In Fig^b-e) we present Husimi distributions for 
several Floquet states. The dimensionless Planck con- 
stant ft = 0.2 is in a range, when it is possible to estab- 
lish a correspondence between different Floquet states 
and the invariant manifolds of the mixed phase space for 
the classical limit. Each plot carries a symbol, which 
shows the corresponding location of the quasienergy of 
that state in FigEt. The states (b-d) are located in vari- 
ous regular phase space regions, while state (e) is located 
inside the chaotic layer. 

The current for an initial condition \tp(to)) depends in 
general on the initial phase to- Note that in the classical 
case such an initial condition may also lead to some de- 
pendence of the classical current on to, since the initial 
distribution may overlap with different regular transport- 
ing manifolds. However, if we start with a cloud of parti- 
cles exactly located inside the chaotic layer, the asymp- 
totic current will be independent of to for any choice of 
the distribution function over the chaotic manifold. This 
is not true for the quantum case where the current may 
even change its sign with the variation of to- It is a 
consequence of the linear character of the Schrodinger 
equation |ljfl. We will first discuss the results obtained 
after averaging over the initial phase to. Then we can 
assign a unique current value, J = 1/T J Q T J(to)dto, for 
fixed parameters of the ac-field, E±, E 2 , and 9. 
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FIG. 3: (a) The average current J (in units of the recoil 
momentum) vs 6 for different amplitude values of the second 
harmonic, E2: 0.95 (pointed line), 1 (dashed line) and 1.2 
(solid line). Insets: relevant details of the quasienergy spec- 
trum versus 8 in the resonance region for E2 = 1 (top right) 
and E2 = 1.2 (bottom left). The parameters are E\ — 3.26 
and lo = 3. (b) The average current J (in units of the recoil 
momentum) vs 9 for E\ — 3, E2 = 1.5 and w = 1. 



Figure a) shows the dependence of the average cur- 
rent on the asymmetry parameter 9 for the initial condi- 
tion \ip) — |0). The average current J shows the expected 
symmetry properties J(9) = — J (9 + ir) = —J (—6). On 
top of the smooth curves we find several resonant peaks 
for E2 = 0.95 where the current value changes drastically. 
Comparing with the quasienergy spectrum, these reso- 
nances can be unambiguously associated with avoided 
crossings between two Floquet eigenstates. The Husimi 
distributions show that one state locates in the chaotic 
layer, and another one in a transporting island. Off res- 
onance the initial state mainly overlaps with the chaotic 
state, which yields some nonzero current. In resonance 
Floquet states mix, and thus the new eigenstates con- 
tain contributions both from the original chaotic state as 
well as from the regular transporting island state. The 
Husimi distribution of the mixed state is shown in Fig^ 7 , 
the strong asymmetry is clearly observed. The regular 
island state has a much larger current contribution, re- 
sulting in a strong enhancement of the current. 

^From an experimental point of view a too narrow res- 
onance may become undetectable due to resolution lim- 
itations. We thus studied how to vary the width of the 
resonance without much affecting its amplitude. It turns 
out to be possible by tuning another control parameter, 
e.g. the amplitude E 2 - We increase this field amplitude 
in order to disentangle the two Floquet states and remove 
the avoided crossing. That will happen for some value of 
E% at 9 = ±7r/2. The details of the quasienergy spec- 
trum around that critical point are shown in the insets 
in Fig[2t- The two quasienergy spectra disentangle for 
E2 = 1 but stay close over a sufficient broader range of 9 
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FIG. 4: Current dependence on the intial phase to and 9. The 
parameters are the same as in Fig^. 

values. Thus the resonances become broader, as seen in 
FigEK- Further increase of E% to a value of 1.2 leads to 
a strong separation of the two spectra, and consequently 
to a fast decay of the amplitude of the resonance. 

The above numerical results for uj — 2 show a maxi- 
mum current value in the resonance region of the order of 
0.06 in units of recoil momenta. In order to increase that 
output, we drive the system into strongest resonance by 
choosing to = 1, since that driver frequency matches the 
oscillation frequencies of particles at the bottom of the 
spatial periodic potential. The result is shown in Fig|3|3. 
We again observe a clear broad resonance, but the max- 
imum current value increases by an order of magnitude 
up to 0.5 in units of recoil momenta. 

As already mentioned, without averaging over to, the 
current depends on the initial phase. However the ob- 
served resonance structures are due to resonant inter- 
action between Floquet states, or avoided crossings of 
quasienergies. These resonances are independent of the 
initial phase to- Indeed, in Fig0]we plot the nonaveraged 
current as a function of both 9 and to- While the smooth 
background is barely resolvable with the naked eye, the 
resonances are clearly seen, and their position is not de- 
pending on to, while their amplitude does. That implies 
that one can further maximize the resonant current by 
choosing proper initial phases to, reaching values above 
the recoil momentum. 

Note that our approach is very different from a re- 
cently proposed modified kicked rotor model 0, where 
a biased acceleration appears for an initial condition with 
preassigned nonzero velocity. Here we study the regime 
of stationary current for the general case of an initial 
state of zero momentum for a model which has a finite 
current also in its classical counterpart [rjH^. 

In summary, we have studied the mechanisms of aver- 
age current appearance in driven quantum systems with 



broken symmetries. The key source of such directed 
transport is the desymmetrization of Floquet states. A 
peculiarity of the quantum ratchet is the strong depen- 
dence of the current on the initial phase of the applied 
field. Moreover, we found quantum resonances induced 
by avoided crossings between Floquet states which en- 
hance the current drastically. Optimizing the drive fre- 
quency, amplitude and initial phase, resonant currents 
easily reach the recoil momentum value and should be 
experimentally observable using driven cold atoms in op- 
tical lattices. 

We acknowledge fruitful discussions with F. Renzoni 
and M. Weitz. 
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